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Abstract 

The density functional theory originally developed by Hohenberg, Kohn and 
Sham provides a rigorous conceptual framework for dealing with inhomoge- 
neous interacting Fermi systems. We extend this approach to deal with in- 
homogeneous interacting Bose-condensed systems, limiting this presentation 
to setting up the formalism to deal with ground state (T = 0) properties. 
The key new feature is that one must deal with energy functional of both 
the local density n(r) and the local complex macroscopic wavefunction $(r) 
associated with the Bose broken-symmetry (the local condensate density is 
n c (r) = |<l>(r)| 2 ). Implementing the Kohn-Sham scheme, we reduce the prob- 
lem to a gas of weakly-interacting Bosons moving in self-consistent diagonal 
and off-diagonal one-body potentials. Our formalism should provide the basis 
for studies of the surface properties of liquid 4 He as well as the properties of 
Bose-condensed atomic gases trapped in external potentials. 
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I. INTRODUCTION 

This article develops a density functional formalism for dealing with superfluid 4 He, 
with the ultimate goal of understanding the role of Bose-condensation in the low density 
surface region. The use of neutron scattering to study the atomic structure and dynamics 
of superfluid 4 He has been an ongoing research effort at Chalk River since the early fifties. 
This world famous programme grew out of and was nurtured by the atmosphere which Bert 
Brockhouse helped to create. On a more personal note, I have always enjoyed being invited 
to give seminars at McMaster, partly for the stimulating discussions with research colleagues 
but also because I knew that Bert would be in the audience and after it was over, he would 
come over and give some encouraging comments. It is with great pleasure that I dedicate 
this article to Prof. Bert Brockhouse. 

At a phenomenological level, the low temperature properties of bulk superfluid 4 He are 
well understood in terms of Landau's picture of a weakly-interacting quasiparticle gas of 
phonons and rotons. At a more microscopic (atomistic) level, current discussions of super- 
fluid 4 He can be divided into two broad classes. One of these is the field-theoretic analysis, 
which is built on the fundamental role of Bose-broken symmetry. The superfluid phase co- 
incides with the appearance of a macroscopic wavefunction given by the finite expectation 
value of the field operator $(r) = (^(r)) This approach goes back to the pioneering 
work of Bogoliubov and Penrose ||, but was first formulated in a systematic way by 
Beliaev [Q and Hugenholtz and Pines ||. The current status of theories of superfluid 4 He 
based on the key role of the Bose order parameter $(r) is reviewed in a recent book ||. 

A second kind of microscopic theory is built on the use of variational many-body wave- 
functions for the ground state and the low-energy excited states of liquid 4 He. Such many- 
body wavefunctions were first introduced by Bijl |7j and later by Feynman ||, and in their 
current form have become very sophisticated (the correlated basis function approach ||). 
However, while such ground state many-particle wavefunctions (of the Jastrow-Feenberg 
type, for example) lead to very good estimates of the condensate fraction (about 10% at 
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T = 0), the role of a Bose order parameter is not exhibited very explicitly (see chapter 9 of 
1). 

The motivation of the present paper is to set up a formalism which can deal with the sur- 
face properties of superfluid 4 He, taking Bose-condensation into account from the beginning. 
Formally this means we need a theory of spatially inhomogeneous Bose-condensed system. 
In the last decade, there has developed a considerable literature on superfluid 4 He with 
free surfaces (including droplets and films) based on generalizations of the above mentioned 
correlated basis function approach [H| . However, as with the case of bulk liquid 4 He, such 
treatments of surfaces give little (if any) emphasis to the role of Bose-broken symmetry or 
Bose-Einstein condensation. On the other hand, the many-body wavefunctions which have 
been developed [[Uj to describe the free surface of liquid 4 He at T = appear to give realistic 
estimates of the density profile n(r) in the surface region. From these, one finds accurate 
values of the binding energy of 3 He and spin-polarized H atoms bound to the surface of 



superfluid He [12 



Such correlated basis function approaches involve a very heavy computational effort when 
dealing with inhomogeneous systems. An alternative, much simpler theory has developed 
for free surfaces of liquid 4 He based on a density functional theory ]nj. This approach is 
loosely inspired by the very successful theory of inhomogeneous interacting electron systems 
and it is useful to recall some aspects of the latter theory. Almost 30 years ago, Hohenberg 
and Kohn (HK) gave a rigorous formulation |14j] of the ground state properties of interacting 
Fermi systems, in terms of a functional of the local density n(r) = (^(r)i/j(r)). Kohn and 
Sham |15j (KS) used the two exact HK theorems to implement the HK formalism in terms of 
finding the single-particle energies and eigenstates of free Fermions moving in an appropri- 
ately defined self-consistent fields. This Hohenberg-Kohn-Sham (HKS) density functional 
formalism is now the accepted way of dealing with inhomogeneous Fermi systems, build- 
ing on our extensive knowledge of the ground state properties of homogeneous interacting 
Fermi systems. The HKS approach has been generalized to deal with normal systems at 
finite temperatures |L6| as well as BCS superconductors |17| . In the latter case, one works 
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with functionals of the local density n(r) and the local anomalous (or off-diagonal) density 
A(r) = (■0|(r)-?/'|(r)) describing spin-singlet Cooper pairs. Finally, a density functional for- 
malism has been developed for dealing with time-dependent quantities (excited-states and 



linear response functions) of both normal [^] and superconducting Fermi systems ||T9 . 

As we have noted, recent density functional theories of superfluid 4 He with surfaces make 
analogies to the above HKS theory of inhomogeneous Fermi systems. However, there has 
apparently never been a careful study of how to use the HKS ideas to give a rigorous basis 
to a theory of inhomogeneous Bose-condensed liquids, analogous to what has been done for 
inhomogeneous BCS superconductors |17]]. In particular, all current T = density functional 
theories of superfluid 4 He 0] simply assume that the energy functional only depends on 
the local density n(r). There is never any reference to the possibly equally important role 
that the local condensate density n c (r) might play, or more generally, the local macroscopic 



wavefunction $(r) = (V>(r)) = J n c (r)e lS<yr \ Needless to say, it is the finite value of the 
order parameter $(r) which characterizes the superfluid phase below T\ = 1.YJK. To avoid 
confusion, we note that the energy in density functional theories of liquid Helium (see third 
paper of ref. |I3|) is often taken to be a functional of the variable \I/(r) = Jti(r)e lS ^, 



which involves the total local density. While this variable is sometimes referred to as a 
"macroscopic wavefunction", it is clearly unrelated to the Bose order parameter $(r) we 
have introduced above. The insufficiency of theories based on functionals of only n(r) has 



been made especially obvious in recent work [2(| which points out that the low density 
surface region of liquid 4 He corresponds to a dilute inhomogeneous Bose gas with 100% 
Bose condensation (at T = 0). In this surface region, the key function is $(r), with the local 
density being determined by it, namely n(r) = n c (r) = |$(r)| 2 . 

In the present paper, we formulate a density functional theory of the HKS kind for the 
ground state properties of an inhomogeneous interacting Bose-condensed fluid. The key 
new element in our analysis (following the analogous case of BCS superconductors [[H]]) is 
the realization that one must work with functionals of both n(r) and the (complex) order 
parameter $(r). The key theorems of HK and the methods of proof are, of course, valid for 
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Bose as well as Fermi statistics. For brevity, we shall only sketch these arguments when they 
involve the identical steps as in the density functional treatment of BCS superconductors. 
The present analysis of ground state properties can be extended to finite temperatures (free 
energies) following the analogous discussion for BCS superconductors [[H]]. This will be 
reported elsewhere. 

What the present paper accomplishes is to give a formally exact scheme for dealing 
with inhomogeneous Bose-condensed fluids which should ultimately provide a platform for 
specific calculations. In applying the present formalism, one must introduce approximations 
for the correlation energy functionals. This is a separate question, with specific problems 
associated with the anomalous long-range correlations in Bose-condensed systems f2l| , and 
is not treated here. 

In this paper, we mainly use the surface region of superfluid 4 He as an example of an 
inhomogeneous Bose-condensed system. Very recently, Bose-condensation has been finally 
achieved |22|] in a dilute gas of 87 Rb atoms below 200 nK, using laser and evaporative cooling. 
This gas was trapped in a harmonic potential well and as a result, both n(r) and n c (r) are 
highly inhomogeneous. Our present formalism gives a natural basis for generalizing the 



currently available Hartree-Fock-Bogoliubov approximations for such inhomogeneous 
weakly interacting Bose-condensed gases. 

II. HOHENBERG-KOHN FORMALISM 



Our starting Hamiltonian is defined as (compare with [17 

V 2 



H VtV = J drft(r 



2m 



V»(r) 

H — J dr J dr'ffi(r)ffi(r')v2(r — r')^(r')^(r) 
+ J (fru(r)-?/^(r)^(r) 

+ J dr[r](r)ft(r)+r]*{r)4)(r)] (1) 
Ho-tiN + Vz + Vt + VsB. (2) 
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Throughout the analysis, the two-particle interaction w 2 (r — r') is assumed to be fixed. Thus 
H v ^ in (|l]) depends only on: (a) the choice of the external diagonal single-particle potential 
v(r), which couples to the density; and (b) the external off-diagonal (or symmetry-breaking) 
potential rj(r), which couples to the field operators ^(r) and ?/^(r). The first step in the HK 
approach is to note that the ground state |^) of H VjV is a functional of these external fields 
v(r) and n(r), and hence so are the following groundstate expectation values: 

n(r) = (^|^ f (r)^(r)|*) 
$(r) ee (#|$(r)|#) 

$*(r) = (tt|$t(r)|*). (3) 

We note that because of the symmetry-breaking field Vsb i n (UK the ground state eigenstate 
of H VtTI allows $(r) and $*(r) to be finite. Introducing a symmetry-breaking (off- 
diagonal) perturbation as in ([]]) is the standard method [|1],(4] of dealing with the appearance 
of Bose-condensation in an interacting system. It involves treating the condensate as a 
"reservoir" of atoms and thus leads to number non-conservation and to a groundstate where 
$(r) can be finite. The underlying physics of this broken symmetry is the same for both 
homogeneous and inhomogeneous systems and is discussed most clearly in Ref. Q. We 
recall that the essential physics for dealing with BCS superconductors |L7| also involves 
such symmetry-breaking number non-conserving states. 

The second step of HK is to note that, using their famous reductio ad absurdum argument, 
that up to additive constants, one can prove 

v(r) and n(r) are unique functionals of n(r) and $(r). (4) 

Since v and n fix H VtJ] and thus also we can finally conclude that 

l^) is a unique functional of n(r) and 3>(r). (5) 

In turn, it follows from (|5|) that the expectation value (^l-ffo — uN + V^l^) is a universal 
functional of n(r) and $(r). That is to say, there is no explicit dependence of on the 
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specific forms assumed for v(r) and ry(r), since these can be expressed as functionals of n(r) 
and $(r), as stated in (f|). Summarizing this train of argument, we conclude that 

F[n(r),^(r)] = (^\H -fiN+V 2 \^) (6) 

is a universal functional of n and $, valid for any number of particles and any external 
potentials v(r) and r)(r). As in the HK formalism for Fermi systems, the universal functional 
F[n, $] will play a central role in our subsequent analysis. A key problem, of course, will 
be to find some appropriate approximation to this universal functional F[n, $] in superfluid 
4 He. 

Following HK, it is useful to define, for given v and 77 potentials, the energy functional 

E VjV [n, $] = F[n, $] + J drv(r)n(r) 

+ J rfr[77(r)$*(r) + r?*(r)$(r)], (7) 

where the densities n(r), $(r) and $*(r) are defined in ([$]). For simplicity of notation, we 
shall generally show functionals as depending only on $, but in fact, they depend on both 
$ and $*. Following HK, we have a variational principle, i.e., one can prove that E VyT) [n, $] 
is a minimum at the correct values of the densities n(r) and $(r) produced by the external 
potentials v(r) and rj(r), i.e., 

8E v , v [n,$] = Q 5E ViV [n,$] = Q 
8n(r) 1 ' 5<&(r) 



III. KOHN-SHAM PROCEDURE 



Following Kohn and Sham [13], there is a clever way of finding the correct values of 
n(r) and $(r) (which, according to (^), minimize S^fn, $]) by solving a simpler auxiliary 
problem for which the HK theorems in Section II are also valid. To understand the logic of 
the KS procedure in the context of our present problem, let us consider an auxiliary system 
Hamiltonian defined by 



2m ^ 



+V S $\$\+ J drv s {v)^{v)^{r) 

+ J dv[r] s {v)^{v) + ry*(r)^(r)], (9) 

where the interaction V, is a part of Vi in (|1]) and (|2]), to be specified later. All the HK 
results of Section II apply to flj|). In particular, if we denote the ground state of £T* a as 
l^s), then the densities 

n s (r) = (*,|^+(r)^(r)|* a > 

$ a (r) = <*.$(r)|tf a ) (10) 

are unique functionals of the external fields v s (r) and r) s (r). In turn, v s and r) s and hence 
\^f g ) can be shown to be unique functionals of n s (r) and 3>«(r). Thus we conclude that 

F.[n a (r), $ s (r)] = (^ s |# - /xiV + V a \V 9 ) (11) 

is a universal functional of n s (r) and $ s (r). Finally, we can define an energy functional of 
this auxiliary system 

K., v .[n,$] = F s [n,$] + J drv s (r)n(r) + J drfo.(r)$*(r) + ^(r)$(r)], (12) 

which will be minimized by the correct values, n(r) = n s (r) and $(r) = $ s (r), for this 
system. 

The whole point of introducing this auxiliary system defined by (0) is that: 

(a) It will be easier to solve than the actual system described by ([[]). 

(b) By a judicious choice of the external fields v s (r) and r) s (r), the densities n s (r) and $ s (r) 
of this auxiliary problem can be made to be identical to those of the real system. Since 
F[n(r),$(r)] in ([]) is a universal functional of only n(r) and $(r), this means that 
we can evaluate it using results for n(r) and $(r) obtained from solving the auxiliary 
system. 
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In applying the KS procedure to superconductors |17j , one uses a non- interacting gas of 
Fermions moving in external one-body and pair potentials as the auxiliary model system. 
In our interacting Bose system, this would correspond to setting V^f^,^] i n ® to zero. 
The problem with this choice is that for a non-interacting Bose gas moving in given external 
potentials v s and r] s , and at T — 0, one has complete Bose-Einstein condensation (see 
discussion after (|35| ) for more details). This implies that n(r) and n c (r) = |$(r)| 2 are equal, 
even though we know that in any interacting Bose system, the local condensate density n c (r) 
is less than the local total density n(r). This problem is not addressed in density functional 
theories [|l3j of superfluid 4 He based on functionals of only the density n(r). We recall that 
such theories usually start with the kinetic energy of an inhomogeneous non-interacting Bose 
gas with a density profile n(r) identical to the fully-interacting system. Such a kinetic energy 
functional implies that n c (r) = n(r), which would not appear to be a very good starting 
point for describing superfluid 4 He. 

In order to define our auxiliary Bose system in (||), we first introduce the usual decom- 
position of Bose quantum field operators || 



-0(r) = $(r) + ^(r) 



(13) 



where $(r) is defined in (|3]). The non-condensate field operators ip(r) and ^(r) satisfy Bose 
commutation relations. Using (|13"D , the two-particle interaction V 2 in (fj) can be rewritten 

as 



'\|2 



V 2 = -JdrJ dr'v 2 (r - r')|$(r)| 2 |$(r 

+ JdrJ dr'v 2 (r-rO|$(r)| 2 $(rO^V) 
+ JdrJ rfr'^; 2 (r-r / )|<I>(r)| 2 <I)*(r , )^(r , ) 

+ / dr f dr'v 2 {r-r') [^(^^(^^(r')! 2 + ^ f (r)^(r')$*(r')$(r 



+ ^ t (r)^ t (r / )$(r')$(r) + ^(r)^(r')$*(r')$*(r) 
+ f dr [ dr'v 2 (r-r') [^(r')^^')^^)**^) + ^(r')^^)^^')*^) 
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+ \l dr I dr ' V2 ( r ~ O^WVM 1 "')^^)- ( 14 ) 

If all the atoms were Bose-condensed, only the first term in ([L4]) would be important. If the 
system is not Bose-condensed, then only the last term in ( |14|) is present. In the well-known 



Bogoliubov approximation [fM|[25| for a dilute, weakly interacting gas, in which almost all 
the atoms are Bose-condensed, one only keeps terms up to quadratic in the non-condensate 
field operators if) and i[)> (since it is assumed that, in an average sense, if) -C 3>). That is 



to say, the last two terms in (HP are higher order and hence omitted. A feature of this 
Bogoliubov approximation is that the resulting Hamiltonian can be diagonalized exactly 
(see below). 

We now define what we shall call (for want of a better term) the exchange-correlation 
energy functional F xc [n, $] by writing @ in the form 

F[n, $] = F s [n, $] + ^JdrJ dr'v 2 (r - r')n(r)n(r') + F xc [n, $], (15) 

where F s [n, $] is the energy functional of the auxiliary system defined by ( |TTD with the 
interaction 

V s = J drj dr'v 2 (r - r / )|$(r)| 2 [$(r , )^ t (r') + $*(r')^(r')] 

+ \l dr I dr W r_ 0[2$V) $ ( r )^ r )<K r/ ) 

+$(r')$(r)^ t (r)^ t (r') + $*(r , )$*(r)^(r)^(r')], (16) 
and subject to potentials v s (r) and r] s (r) [see (|TJ)] chosen such that the density n s (r) and 



order parameter <& s (r) are identical to those of the full system. As usual ||l^Jl5H , it is useful 
to separate out the total Hartree energy contribution as we have done in (15). Writing this 
contribution out more explicitly for a Bose-condensed system, we have 

(V H ) = \j dr J rfr '^(r - r')[|$(r)| 2 |$(r')| 2 + 2|$(r)| 2 n(r') + n(r)n(r% (17) 

where the non-condensate local density is defined by 

n(r) = (ft(r)ip(r)) = n(r) - |$(r)| 2 . (18) 
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The last two terms in (|17D come from the ip^tp and (ip^ifj) 2 terms in (|H|) . 

Calculating the variational derivatives in (§) using $] in (|i~5|) , one finds 

ggjryg] | 5F xc [n,$] 
5$(r) 5$(r) 

where the Hartree field is defined as 



7y(r) = 0, (19) 



v H {r) = J dr'v 2 (r - r')n(r'). (20) 



Similarly, using (Q) for the auxiliary system defined above, one finds 

5F a [n,$] 



5n(r) 
SF s [n,$] 



+ v s (r) = 

+ Vs(*) = 0. (21) 



<J$(r) 

Combining the results in Cp~9|) and (|2l|) , we conclude that the density and order parameters 



of the auxiliary system will be identical with the actual system if 

v s (t) = v(r) + v H {r) + ^fe^ = v a [n, $] 

%(r) = r ? (r) + ^^l = r ?s [n,$]. (22) 

This gives u s and ?7 S as explicit functionals of n(r) and $(r), once we have decided on a 
specific form for the functional F xc [n, $]. 



IV. BOGOLIUBOV GAS AS AUXILIARY SYSTEM 

We now turn to the auxiliary system defined by (§), with V s given by dl6|) . Using (|13|). 
one finds that 

= / dr**(r)[Jte(r) +17.(1:)] + / ^ s *(r)$(r) 
+ y dr^ f (r)[L<l>(r) + ry s (r)] 
+ 1 dr[$*(r)L + 77 s *(r)]^(r) 
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+ J dr^ f (r)L^(r) 



+-J drj drV 2 (r-r') [2$*(r / )$(r)^ t (r)^(r / ) 
+$(r / )$(r)^ t (r)^ t (r / ) + $*(r')$*(r)^(r)^(r') 



where the operator L is defined by 



V 2 \ 
"2^ + «.(')" A* • 



(23) 



(24) 



This Hamiltonian is similar in structure to the one one obtains in an inhomogeneous weakly 
interacting Bose gas at T — 0, which has been extensively studied in the literature p4] , |25| . 
We solve it using similar techniques. In order to diagonalize (0), we first eliminate the 
terms linear in ip and i\y by requiring that $(r) satisfy the equation 



L$(r) +r} s (r) = 0. 



(25) 



Eq. ( |25| ) is a sort of generalized Gross-Pitaevskii equation for but now in the context of 
density functional theory rather than for a dilute Bose gas. This connection is easily seen 
by setting F xc [n, $] in ( Jl5| ) to zero, in which case ( [22] ) simplifies to 



v s (r) = v(r) + v H (r) 
r) s {r) = rj(r), 



(26) 



and (|25|) reduces to 



V 2 r 
— — + v(r) - /i + J dr'v 2 (r - r )n(r ) 



$(r) +?7(r) = 0. 



(27) 



Setting the external fields v and i] to zero, we recover the well-known Gross-Pitaevskii 



equation p6| , p4]| for a dilute inhomogeneous gas at T = 0. In that case, since all the atoms 
are Bose-condensed, the density n(r) in (p7|) can be approximated by n c (r) = |$(r)| 2 , and 
then (|27| ) is a closed non-linear Schrodinger equation (NLSE) for $(r). 

Assuming that $(r) satisfies (123), our auxiliary system Hamiltonian (P3"D reduces to 
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H 



Vs,v s 



dr?7*(r)$(r) 



+ J drfo(r) 



v 2 



+ dr dr'v 2 {r - r')<f>*(r')<$>(r)ft(r)4>(r') 



+ -j dr j dr'v 2 {r-r') [$(r / )$(r)$ + (r)$ t (r / ) + $V)$*( r M r M r ') 



(28) 



The Hartree contribution is contained in u s (r) [see (p2|)1. The third line in (B3) gives the 
exchange term, while the fourth line involves the anomalous contributions involving i/j'tp' 
and ipip characteristic of a Bose-condensed system. The quadratic expression given by 
can be diagonalized by the usual Bogoliubov transformation [EH] 



$(r) = l u j( r ) a j - v*(r)a]} 



(29) 



where the new "quasiparticle" operators a 3 - and ft] satisfy Bose commutation relations. 
One finds the amplitudes Uj(r) and Vj(r) are given by the generalized Bogoliubov coupled 
equations: 



-2^ + Vs[v) -» 



Uj(r) + J dr'v 2 (r - r')$(r)$*(r')uj(r') 

- / dr'v 2 (r - r')$(r)$(r> j (r') = £>j(r) 
Uj-(r) + J dr'v 2 (r - r / )$*(r)$(r')wj(r / ) 

-J dr'v 2 (r - r')<5>*(r)<5>*(r') Uj (r') = -E jVj ( 



(30) 



One can prove that the eigenvalues Ej are real and that one must choose solutions such that 
Ej > 0. For details of how (|30"D is derived, we refer to a similar calculation by Fetter [25 



for the case of a contact interaction v 2 (r) = vo5(r). However, we note that long-range tail 
of the He-He interatomic potential is very important when dealing with the surface region 



of liquid Helium |20 |. 

With Uj and Vj given by the solutions of Q50]), the Hamiltonian in (E8| ) can be shown to 
reduce to 
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Hl Vs = J dr V * s (r)<f>(v) - £ Ejf rfr|^(r)| 2 + £ i'. j<\ a r (31) 

which describes a non-interacting gas of quasiparticles of energy Ej. The ground state \^ s ) 
of this Hamiltonian is defined by <%|\I/ S ) = 0. Thus the ground state expectation value of 
(|31~|) is given by 

E^ Vs = J di^(r)*(r) - c/r|^(r)| 2 (32) 

Using (|29|) , the non-condensate local density in (|18D is given by (T = 0) 

n(r) = <tf 8 |$t( r )$( r )|* a ) 

= EK-WI 2 (33) 

i 



Inserting ( pT|) into (]T2[) and using (|25|) and (|i~8D, one finds after a little algebra that 

F s [n,$] = |dr$*(r) - // $(r) - ]T | dr|^-(r)| 2 + i;,(r)] . (34) 

L J j 

We recall that the Hartree contribution was separated out in ( |T5| ) and consequently it is 
not contained in /f* defined in Thus it is not included in the energy eigenvalues 

£y given by fl30|), but rather appears as a separate contribution from the diagonal potential 
v.(r) in (H. 

The key feature of the above results is that one can have a depletion of the condensate, 
as shown by the finite value of n(r) in ([33]). Thus the auxiliary system defined by (|23|) and 
( p2| ) can be used to find both $(r) and n(r) [using (p5|) and (0)] even when n c (r) = |$(r)| 2 
and n(r) are quite different (as in superfluid 4 He). 

If we had chosen a non- interacting Bose gas as our KS reference system [i.e., set V s = in 
(§)], the last term in ( p3|) would be absent. The linear terms in ip and ^ can be eliminated 
as before by requiring that $(r) satisfy (|25|). Then H* a is easily diagonalized, with Uj(r) 
given by the solution 
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and Vj(r) = 0. It immediately follows from (^) that such a KS reference system leads 
to no depletion, with n(r) = 0. This means that using a free Bose g reference 
system inevitably leads to n c (r) = |$(r)| 2 and n(r) being identical, no matter what we 
choose for the functional F xc [n, <£>]. This apparent "insufficiency" of the non- interacting 
Bose gas is somewhat surprising and deserves further study. It seems to be associated with 
the well-known fact that a weakly interacting Bose-condensed gas has qualitatively different 
properties than an ideal Bose gas. 

It is convenient at this point to summarize the various steps in the KS procedure: 

(1) One chooses some approximation for the "exchange-correlation" functional F xc [n, $] 
defined in (|15|), giving it as an explicit functional of the local quantities n(r) and $(r). 



This is the big step containing the "physics" , and has not been addressed in the present 
paper. 

(2) The potentials v s (r) and r) s (r) are then computed using (|^) and given as functionals 
of n(r) and $(r). 

(3) Evaluating v s and r] s using an assumed (or trial) value of n(r), the GP-type equation 
25|) is solved for §(r), i.e., 

V 2 



+ v a [n,$]-(jL 
2m 



$(r) +r] s [n,<f>] = 0. (36) 



(4) Evaluating f s [n,$] using the trial n(r) and the solution for $(r) given in step (3), 
the generalized Bogoliubov equations in (|30|) can be solved to determine Ej, Uj(r) and 



Vj {r). 

(5) With these results, one can calculate the non-condensate local density n(r) in (|33D and 
hence finally obtain n(r) from (|I8|) , namely 

^r) = |$(r)| 2 + £Mr)| 2 . (37) 

3 

Using this new expression for n(r), one can go back to step (3) and repeat the proce- 
dure, until self-consistent values of $(r) and n(r) are obtained. 
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(6) These values of n(r) and $(r) are then inserted into the energy functional of the actual 
Bose system of interest, as given by (|7|) and flT5l). 



V. CONCLUSIONS 



As in the case of interacting Fermi systems |To| , |l7j| , it is important to emphasize that while 
the auxiliary Bose system described in Section IV corresponds to a dilute, weakly interacting 
Bose gas, we are only using it to find the local density n(r) and local order parameter $(r) of 
a Bose-condensed liquid. In particular, the quasiparticle excitations which are described by 
the generalized Bogoliubov equations of motion in (^) have, in general, no "direct" physical 
significance in the true system. 

The above procedure gives a well-defined scheme to find (at T = 0) the energy, local den- 
sity and local condensate density in superfluid 4 He with a free surface, once one has chosen 



some explicit approximation for the exchange-correlation functional F xc [n, $] in (|15|). As in 
the case of normal and superconducting metals, a good approximation for this functional is 
the key to obtaining reasonable results using the density functional formalism. We hope to 
discuss this problem elsewhere. In the case of superfluid 4 He, there has been considerable 



work [13| on constructing functionals which only depend on the total local density n(r), in 
which one tries to build in known experimental information (compressibility, ground state 
energy, surface tension, etc.). In developing the equivalent approximations for use in our 
new formalism, the first thing we need to understand better is the ground-state energy of 
bulk liquid 4 He as a function of the density n and the condensate density no- More Monte 
Carlo calculations would be very useful. 

One immediate implication of the present theory is contained in the first term of fl34|), 



which can be rewritten in the form (for clarity, we insert K) 

J dr$*(r) $(r) = J dv^-\V ^f^r)\ 2 + drmn c {r)v 2 s (r) , (38) 



where we have used $(r) = Jn c (r)e lS ^ and mv s (r) = h\/S(r) is the local superfluid 
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velocity (in the ground state, we can set v s (r) =0). In contrast with the first term in ([38] 
currently available density functional theories |13j always start with a term of the form 

n 2 



dr i |V \M r)|2 ' (39) 

involving the total local density n(r). This corresponds to the kinetic energy functional of an 
inhomogeneous non-interacting Bose gas constrained to have the correct local density of the 
Bose liquid under consideration. We believe that the first term in (|38D has a more natural 
as well as more sound theoretical basis when dealing with Bose-condensed fluids, and should 
be used in developing improved density functional theories of superfluid 4 He [^0f . 

At the present time, the only detailed discussion of the properties of a spatially inho- 
mogeneous Bose-condensed system has been for a weakly interacting atomic gas trapped 
in an external potential well p3| . In such systems p^j , the condensate density n c (r) is 



strongly peaked at the center of the trap (due to the macroscopic occupation of the lowest 
quantum state) and has a very different spatial dependence from the non-condensate density 
h(r). Our density functional formalism (extended to finite temperatures) should be useful 
in calculating the properties of such trapped Bose-condensed gases. 
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